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1 INTRODUCTION 



With this paper we continue a series of investigations devoted to multidimensional grav- 
itational models with the aim to obtain exact solutions and to find on their basis some 
observational consequences of extra (hidden) dimensions in our physical world. We also 
study the role of different physical fields in a self-consistent manner. 

Modern theories of field unification assume in general a greater space-time dimension- 
ality than four. Although multidimensional gravity as an approach to field unification 
can be traced back to the famous works of Kaluza and Klein [1,2] in the twenties, to- 
day's increased interest to this field is largely stimulated by studies in supersymmetry and 
some other modern theories [3]; in the field-theoretical limit of such theories gravity is 
described with reasonable accuracy by multidimensional Einstein equations. Studies of 
their solutions can lead to predictions of direct observational manifestations of extra di- 
mensions. Thus, cosmological models predict variations of the gravitational constant G, 
so that observational constraints imply certain limits on model parameters [4,5]. In these 
papers exact solutions as well as relations between the rate of a possible change of G with 
the Hubble parameter, the deceleration parameter and the mean density of the Universe 
were obtained. Another possible window to the multidimensional world is opened by an 
analysis of local effects which could be sensitive to spatial variations of extra-dimension 
parameters [6]. We shall discuss some effects of this sort, in particular, those connected 
with electric charges of isolated bodies [6-8]. Quantum variants of all these models [4] deal 
with problems of quantum wormholes, creation of the Universe via tunnelling, nonsingular 
models etc. 

Here, we consider exact, static, spherically symmetric solutions to the Einstein-Maxwell- 
scalar field equations, with a dilatonic type coupling between the scalar and vector fields, 
in I?-dimensional gravity with a chain of n Ricci-fiat internal spaces and an arbitrary di- 
mension of coordinate spheres (orbits of the group of spatial motions) — see Eq.(6) for the 
space-time structure. In particular, we pay attention to the arbitrariness of d, in the spirit 
of Tangherlini's generalization of the Schwarzschild metric [9], with the hope to shed some 
hght at the role of the 4-dimensional nature of the physical space-time. 

Our approach differs from that applied in some papers on multidimensional black holes 
in that any solutions, not only black-hole ones, are sought. Consequently, the place of black 
holes (if any) in the whole set of solutions, the role of fields for the existence of different 
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types of solutions, as well as the properties of all spherical configurations, become clearer. 
Here, we call a black hole solution a configuration with the central singularity screened 
by an event horizon. In addition, we study the stability of these solutions under small 
perturbations preserving spherical symmetry and arrive at a result of physical interest 
that from the whole set of static solutions only the black-hole ones are stable. 

2 BASIC EQUATIONS AND THE GENERAL 
STATIC SOLUTION 

Let us consider a system with the Lagrangian 

L = + g^'^^M^'N - e^'^F'^^FMN (1) 

for the set of interacting fields arising in the field limit of superstring theories. Here <^ is a 
dilaton scalar field and Fmn = f^M^w ~ Q^Am is an Abclian gauge field interpretable as 
the electromagnetic one. As pointed out in Ref. [10], an electromagnetic field introduced in 
the multidimensional action may seem less aesthetic than a purely gravitational (Einstein) 
action usually considered in Kaluza-Klein theories but it appears that elementary gauge 
fields are necessary also in higher dimensions in order to obtain a realistic grand unification 
theory. 

The corresponding set of field equations is 

V^Vm^^ + Ac^^^F^^Fm^ = 0, (2) 
VNie^^^F^^) = 0, (3) 

RmN — gMNRA./"^ — —TmN (4) 

where Tmn is the energy-momentum tensor 

Tmn = VmVn - ^QmnV^a + e^^'^[-2F^FNA + ^QmnF^^Fab]- (5) 

Capital Latin indices range from to D — 1, the gravitational constant is put equal to 
unity. 

Let us choose a L>-dimensional manifold of the structure 

n 

M = M(^+'') X Ml X ••• X M„; dimM^ = iV^; D^3 + d + Y,Ni, (6) 

i=l 
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with the signature (+,—,...,—), where M^^'^^^ plays the role of the ordinary space-time 
and Mj are Ricci-flat manifolds with the hne elements ds'^^-y i ~ 1, - ■ ■ ,n. We seek solutions 
of the field equations such that M*^^+'^) is a static, generalized spherically symmetric space- 
time with the metric 

dsl^^ = e^^^^^de - e2"(")dti2 - e^^(")dQ2^i (7) 

where c?i^^+i is the line element on a unit {d + l)-dimensional sphere S"^^^, while all the 
scale factors e^» of the internal spaces Mj depend on the radial coordinates i.e., the 
D-metric is 

n 

dsl = gMNdx^dx"" = dsl^^ - Y: e'^'^-^dsly (8) 

1=1 

If we denote 7 = iV_i = 1, (3 = (3o, Nq — 2 and choose the harmonic radial coordinate 

u such that 

n 

a=Y. (5iN,. (9) 

i=-l 

the Ricci tensor components it!j^ can be written in the highly symmetric form {x^ — u) 

n 

i=-l 

= {N>d + 2; /i = 0,---,(i + 2); 

Rl = ... = it:^+^ = (d+l)e-^^-e-^"/3"; (10) 

where the indices aj{bi) refer to the subspace Mj{Mi). 

Let us further assume the dilaton scalar field to be <^ = (p{u) and the electromagnetic 
field to be Coulomb-like: Am — S^Aolu). Then the D-dimensional Maxwell-like equations 
give: 

= qe-^^"^/^ = ge-2^^-^", q = const (charge), 

^ = I det qmn] = exp(2a + 2 N^) = e^'*. (11) 

i=-l 

(the metric determinants of ds(^i),i — 1, . . . ,n, never appear in the equations and may be 
omitted) . 
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Now the scalar and (some linear combinations of) the metric field equations may be 
written in the form 



if" + 2q^X^^-^^'P = 0; (12) 

Np'^ + j" = 0, i^l,...,n; (13) 

(a -/?)"- dV"-'^ - 0; (14) 

^ (15) 

n 

- E - d{d + l)e"'-^f = ip'^ - 2g2e2^-2V (16) 

i=-l 

where Eq.(16) is the (") component of the Einstein equations and forms a first integral of 
the remaining equations (12-15). 

Eqs.(13, 14) are easily integrated to give 

Pi — —^/N + hiU, hi — const, i — 1, . . . ,n; (17) 
e^-" = d ■ s{k, u), k = const (18) 

where for any choice of variables 

sinh ax, a > 0; 
X, a = 0; (19) 

sin ax, a < 0. 

and inessential integration constants have been eliminated by shifting the origin of u and 
rescahng the coordinates in the subspaces M, (their true scales are thus hidden in dsf^^^). 
Certain combinations of Eqs.(12) and (15) are also easily integrated giving 

= Cu/A-2XN+u;, (20) 
7 = {u; + XCu)/A, (21) 

with the function u!{u) defined by 

e^"^ = Qs{h,u + ui); /i,mi = const; Qs{h,ui) = 1 (22) 

where h and C are integration constants; the other constants are defined by 

A^1 + X'^{N + 1)/N, Q^^q'^/N+, N+ ^ {N + 1)/{2AN). (23) 
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Since, by our choice of the origin of u, e^ — > oo for w — > 0, the value u — corresponds to 
spatial infinity. The last condition from (22) is the requirement that 7 = at infinity, i.e., 
dt is a time interval measured by a distant observer's clock. 
The final form of the £)-dimensional metric is 



dsl = e'^dt'' - e 



-27/Ar 



1 




2/d 


du^ 




d ■ s{k, u) 




d'^s'^{k, u) 



^XCu 



l/A 



(24) 



(25) 



Here and henceforth Yj = 
d+1 



Qs{h, u + ui 

J2i=i ■ The integration constants are connected by the relation 



, A;^signA; = 2N+h'^signh + —r + —r + ^ihl- 
d Ad 



(26) 



This general static spherically-symmetric solution has (n + 3) essential integration con- 
stants: the scalar charge C, the electric charge q (or Q), the "charges" of extra dimensions 
hi and one of the constants h or k connected by Eq.(26). For d — 1 (i.e., conventional 
spherical symmetry with 2-dimensional spheres) h along with C and Q can be related to 
the mass m which is defined by comparing the asymptotic — > (r — > 00) of e''' with the 
Schwarzschild metric: 

AGm + AC = (g^ + h^signh) , (27) 
The coordinate u is defined in the domain [0, 00) if /i > 0, Mi > or between zero and some 
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> in other cases. The scale factors e^^ = 1 for m = 0. 



3 PROPERTIES OF THE SOLUTION 
AND SPECIAL CASES 

Let us consider some special cases of the solution. 

(a) The electromagnetic field is eliminated when Q — 0. Then we recover the generahzed 
[7] Tangherhni solution [9] ; the transformation 



(28) 



6 



brings it to the form 

dsi = rdi" - [i^v^-"-^/'^ 

{d+l)/d = Y.Niaf + a^ + {a + b)yd+C^/e (29) 

where b — ^i^i- The (o?+3) dimensional part of (29) coincides with the Schwarzschild 
solution when d — 1 (the physical space-time is 4-dimensional) and hi — 0. 

The scalar field in (29) affects the metric only via the constant C in the relation among 
the constants. 

The space-time (29) has a horizon (at r = 2k) only in the simplest case a — 1 = = 
b = C = 0. Thus in our model an electrically neutral black hole can exist only with 
"frozen" extra dimensions (e'^* — const) and with no scalar field. (The latter result is 
well known in conventional general relativity: a massless, minimally coupled scalar field 
is incompatible with an event horizon [11]). In this sense the generalized Schwarzschild- 
Tangherlini black holes may be called trivial. These are examples of the so called " spon- 
taneous" compactification of extra dimensions. 

(b) when A = 0, we obtain the generalized Rcissncr-Nordstrom-Tangherlini solution for 

linear scalar and electromagnetic fields discussed in Ref.[8]. 

(c) The scalar field is "switched off" when X — C — 0, which yields a special case of item 

(b). 

A new feature implied by a nonzero electric charge as compared with item (a) is that the 
constants k and h can have any sign and the function s{h,u + Ui) can be sinusoidal, which 
yields i^max < oo. Physically that means the appearance of a Reissner-Nordstrom repelling 
singularity at the center of the configuration. 

(d) When extra dimensions are absent, we arrive at a solution for interacting scalar and 

electromagnetic fields in (d+3)-dimensional general relativity; for the conventional 
case 0? = 1 it was obtained in Ref.[12] (s. also [14]). 



i=l 
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(e) For the case of 4-dimensional physical space-time (d=l, i.e., in Eq.(5) the space 
j^s+d ^ ^4 ^ R^xR'-xS^), the corresponding solution was obtained in Refs. [6,7,13]. 



Comparing the solutions with and without extra dimensions, one sees that indeed the 
scale factors Pi of the Ricci-flat spaces Mj behave like additional minimally coupled scalar 
fields. 

For nonzero electric charge, i.e., in the general case, there is no choice of integration 
constants such that the extra dimensions are frozen = const). So, in this case we have 
no solutions with spontaneous compactification. The behaviour of the metric coefficients 
for different combinations of integration constants is rather various: thus, for — > oo 
(provided h > 0) 



AN 



-h + XC ^ 
u, 1' u, Pi 



-h + XC , 



u, (30) 



so they may be finite or infinite. Calculations show that the solution has a naked singularity 
&X u — limax or 14 = oo in all cases except 

hi = -k/N; h = k; C = -Xk{N + 1)/N, (31) 

when the sphere u — oo is a Schwarzschild-like event horizon: at finite radius r = e^ of a 
coordinate sphere the integral / exp(a — 'j)du for the light travel time diverges. So, (31) is 
the black hole solution we are looking for. 

In general, four types of field behaviors may be singled out: 

A. If /i < or (and) ui < 0, which may take place only when there is a nonzero charge 
q, the coordinate u is defined between and a certain finite value Wmax = zero. 
s{h,u + Ui) < oo corresponds to a repelling (since goo — > oo) naked Reissner- 
Nordstrom-like singularity where the electric field tends to infinity while the scalar 
field (f and the extra-dimension scale factors Pi remain finite. 

B. The case h > 0, Ui > 0; r — as m — >■ oo. The value w — > cx) corresponds to 
an attracting (^fQo scalar-field dominated naked central singularity where the 
dilaton field (p and the "scalar fields" Pi (at least some of them) become infinite. 

C. The case h > 0, ui > 0; r — > oo as -u ^ oo. This configuration, sometimes 
called a space pocket ( "Raumtasche" by P.Jordan), possesses an attracting naked 
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singularity located at a sphere of infinite radius hidden beyond a regular "throat", 
i.e., a minimum of the coordinate sphere radius e^ which is reached at some finite 
value of the u coordinate. Again at least some of the "scalar fields" are infinite at the 
singularity. The system geometry would be like that of a wormhole if it were regular 
at u — oo. 

D. The special case (31), which is intermediate between cases B and C, corresponds to a 
configuration called a dilatonic black hole. In General Relativity {D — A^d— 1) such 
a solution was obtained for the first time in Ref.[12] (see also [14]) and was recently 
widely discussed in various aspects (see, e.g. [15] and references therein). We will 
also pay some attention to it. 

In the black-hole case (31) only two independent integration constants remain, say, k 
and Q, and the transformation (28) brings the solution to the form 

(1 - 2k/R)dt^ 



ds^ 



D 



'R'-'^V'" dR^ fR\^l^ 



-{l+p/Rf/^^ 



dd+i I l-2k/R \dj 



(32) 



V = ^^ln(l+p/it:), (33) 



P = jQ^ + k^-k. (34) 



"^^^ +R'dn^ + J2ds^ 



(35) 



It is a generalization of the Reissner- Nordstrom solution. For d — 1 the metric takes the 
form 

2 _ (1 - 2k/R)dt^ >„/u)VAN 
"^'^^ {l+p/Rfl^ [l-2fc/i? 

considered in Rcfs. [16] (in other notation) and [13] and is in tm^n reduced to the genuine 
Reissner-Nordstrom metric if one puts, in addition, D = A (no extra dimension), A = 
(no scalar field). For = 1, analysing the asymptotic i? ^ oo, one naturally introduces 
the so-called Schwarzschild mass m expressible in terms of Q and k (or k is expressed in 
terms of Q and m): now Eq.(34) is supplemented hy p = A{G'm — k). 

In this family of black-hole solutions a nonzero scalar dilaton field (33) exists solely 
due to the interaction (A ^ 0). When A = 0, i.e., the if field becomes minimally coupled, 
a horizon is compatible only with ip =const. This conforms with the well-known "no- 
hair" theorems and the properties of the general-relativistic scalar-vacuum and scalar- 
electrovacuum configurations [11,14]. 
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It should be pointed out that the solution (32-34) is a very special case of the general 
solution (2 integration constants instead of (n + 3)). So there should exist very strong 
arguments for an assertion that in multidimensional gravity an actual collapse of a spher- 
ical body should lead to black hole (BH) formation. Such arguments do arise when we 
investigate the stabihty of the static solutions. 

It can be demonstrated that one of the specific potentially observable effects for multi- 
dimensional systems is the violation of the Coulomb law. Indeed, the electric field strength 
E = (F^^Fio)^/^ reads for our general solution (10) in the physical case d = 1 

E ^ {\q\/r'')e-''~^^'^ (36) 

where r{u) — is (as before) the curvature radius of a coordinate sphere (the curvature 
coordinate) and 

n 

a = J2NA- (37) 

i=l 

As seen from Eq.(35), the deviations from the Coulomb law are both due to extra dimen- 
sions and due to the scalar-electromagnetic interaction. 

In the BH case for o? = 1 from the asymptotic (r oo) of the metric (35) we have: 

^ = ^ {l - ^ [(Gm + ^^'{Gm + k)^] + O (1) } (38) 

One can conclude that the magnitude of Coulomb law violation is of the order of the 
gravitational field strength characterized by the ratio Gm/r and depends also on N and 
A. 



4 STABILITY ANALYSIS 

Now we would like to investigate the stability of our static solutions under perturbations 
preserving spherical symmetry (i.e., monopole modes). 

As verified by experience, although the D-dimcnsional setting of the problem is quite 
convenient for finding the static solutions, it makes a stability investigation very awkward. 
The latter is better carried out in a {d + 3)-dimensional formulation. 

Prom the viewpoint of the physical space-time V^+s the scale factors Pi are additional 
scalar fields. Thus in the Lagrangian (1) the D-curvature and the metric are to be written 
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explicitly in terms of (3i, after which the Lagrangian acquires the characteristic form for a 
scalar-tensor theory of gravity: 

L = e^ _ af^a^ + J2 fk,M + V'Vm - e'^'^F'^^F^,] (39) 

where corresponds to the metric g^i,, a is defined in (37) and it is assumed that 

<p = cpix"), F^, = F^.(a;°); Fmn = for M or at > d + 2. 

The factor e*^ emerges from the £)-dimensional determinant. 

It is helpful to pass from g^^, to the conformal metric g^,, defined by 

g,^^e-'^'^'^'^g,. (40) 

The Lagrangian acquires the form 

L = R^'^ + ^a'V,, + ^ 7V,A,mA''^ + - e'^/^'+'^^F'^PF^f,, (41) 

where we have denoted d' — d+1, R^'''^^^ corresponds to g^^, and the indices are raised and 
lowered using g^^. 

Taking g^^, in the standard static, spherically symmetric form 

d4+3 = g^^udx^dx" = e^^dt^ - e^'^du^ - e^'^dfl^^^ (42) 

and choosing the harmonic coordinate u, i.e., putting a — {d + 1)P + one rather easily 
restores the general static solution (20)-(24). It is noteworthy that the coordinate u is 
harmonic with respect to both the D-dimensional metric gMN and the (d + 3)-dimensional 
metric g^^, i.e., n^iw = ^d+su — 0, but not with respect to the {d + 3)-dimensional part 
gfj,i, of the D-metric. 

The exphcit form of the metric (42) is easily obtained from (24), (25) using the trans- 
formation (40). 

Now let us investigate small deviations from the static configurations 

5ip{u,t), S(3i{u,t), 5gi^^{u,t), 5F^^{u,t), (43) 

preserving spherical symmetry, i.e. monopole modes. Thus dynamical degrees of freedom 
are restricted to the scalar field and the scale factors which in the 4-dimensional repre- 
sentation behave as effective scalar fields. For simplicity let us assume that there is only 
one internal space: 

D -d-^^ N -d>Q] Ns^ ■■■ ^0. (44) 
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In what follows, with no risk of confusion, we will omit the overbars at the symbols 
q;,/3,7. 

The next step is to choose the frame of reference and the coordinates in the perturbed 
space-time. Evidently this choice may be carried out by prescribing certain relations among 
the perturbations. Following [17], we would like to choose the so-called central frame of 
reference, where coordinate spheres of fixed radii are at rest, and the radial coordinate 
is taken such that the numerical values of these radii are the same as those in the static 
background configuration with the metric (42). Thus we postulate 6f3 = 0, or, as it is 
sometimes more convenient to use the radius r = (coinciding with the Schwarzschild 
radial coordinate), Sr = 0. 

The perturbed metric functions j{u,t) and a{u,t) are taken in the form 

'y{u,t) — ^{u) + 5^{u,t); a{u,t) — a{u) + 5a{u,t), (45) 

and similarly for Lp{u,t), Pi{u,t) = fl{u,t) and F^^iu^t). The perturbed Maxwell-like field 
is defined by AQ{u,t). 

Integrating Eq.(3), we get 

where q is unperturbed since we study only dynamical perturbations rather than changes 
of integration constants and r is unperturbed by the choice of the frame of reference. Prom 
the metric field equations we obtain equations for S/j, and 5(p: 

r^'^'Sil-Sfj," - fj,'{Si -Sa') + 2fj,"{Sa-w) =0, (47) 
r'"^'S(fi-5ip" -cp'iSy -Sa') + 2ip"{Sa-w) =0, (48) 

w = XS(p+{N- d)6i^/d' (49) 

where /i', /i", </?' and (f" are static functions from the background solution and 5^' and 5a 
are defined in (45). The Q^component of the metric field equation is easily integrated in 
t: 

d'p'5a = ~ y ^ V ''^/^ + V'S<P + F{u), (50) 
and the difference of the and components gives: 

d' ofr. f . A {N-d){N+l) ,^ 

-P'{Sa' + 57') = ^ 2> -l^^l^ + (51) 
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Taking 5a and 5^' from (50) and (51) and substituting them into (47) and (48), we arrive 
at coupled wave equations for 5// and 5ip: 



d' 



{N -d){N+l) (rip' fx'' 



d' 



d' 



= 2ip"w] (52) 
= 2fi"w. (53) 



Our static system is unstable if there exist physically allowed solutions to Eqs.(52) and 
(53) growing at i — > oo. A separate problem is to define which solutions should be accepted 
as physically allowed ones. Let us join the approach of Ref.[17] and require 



5// — > 0, 5ip ^ Q for u 







at spatial infinity r — > cxd and 



5fi/ fi \< oo, I 5ip/ip |< oo 



(54) 



(55) 



at singularities and horizons. These requirements provide the validity of linear perturbation 
theory over the whole space-time, including the neibourhood of the singularities. We 
must in addition forbid energy fluxes to our system from outside; however [17], such a 
requirement constrains only the signs of the integration constants and does not affect any 
conclusions. 

There are some cases when the set of equations (52), (53) simplifies, i.e., decouples or 
reduces to a single equation with one unknown function: 

1. the dilaton field is absent: A = 0, = = 0; 

2. there are no extra dimensions: // = (5/x = 0,A/^ = l; 

3. some combinations of (52) and (53) decouple. 

Due to [18], the third possibility is realized for black-hole solutions when d = 1\ how- 
ever, this is probably not the case for arbitrary d. Here we confine the study to the first 
variant when the system dynamics is driven by extra dimensions. We shall verify that the 
arbitrariness of d does not affect the stability conclusions. 

Separating variables in (53) and transforming fi and u to obtain the normal Liouville 
form of the mode equation. 
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5// = e'%(x)//'/^ X ^ - j r'^' {u)du, (56) 
we obtain the Schrodinger-like equation 

y,, + - V{x)]y = (57) 

with the effective potential 



V{x) ^ 



2d' 



dT^ V 7^ j ^ [ r(''+3)/2 j ~~d^^ 



(58) 



where (i' = + 1, again. 

Our static system is unstable if there exist physically allowed solutions of (57) with 
fl^ < (negative energy levels in the language of the Schrodinger equation). 

V{x) is a rather complicated function of x, undetermined explicitly, so that Eq.(57) 
can hardly be solved exactly. However, the main conclusions on the system stability can 
be drawn on the basis of the qualitative behavior of V{x) and the asymptotic forms of the 
solutions to Eq.(57). The latter can be classified according to different behaviors of the 
background static solutions for different values of the integration constants (the notations 
A-D coincide with those in Section 3) . 

A. /i < or (and) Ui < 0, the case of a repelling Reissner-Nordstrom singularity. In 
the limit u — >• Wmax [18] 



X ~ (u^.. - uf+'Z'' , (59) 

so that V has a negative pole. The boundary conditions (54), (55) take the form 

y — > at X — > oo, (61) 

^-(iv+i)/(4iv+2)|^|^^ at x^O. (62) 

At spatial infinity {x — > oo) V{x) behaves hke const/x''+^ (in this and all the other cases). 

We see that both the potential and the boundary condition near x — are independent 
of d and coincide with those derived in Refs. [6,18]. Therefore the conclusion is also the 
same. It looks as follows. 
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The asymptotic forms of the solutions to Eq.(57) for < at x — > and x — > oo are 

x^oo: y = Cie-I^l^ + Cae'^l^, (63) 

x^O: y = C3X^l+C4X^^ si,2 = |(l±2^). (64) 

The condition (63) is satisfied for any C3 and C4. Therefore the solution y{x) chosen on 
the basis of (61), i.e., that with C2 = 0, is admissible and realizes the instability of our 
static system. 

In other words, in the potential well corresponding to the bare singularity there exist 
arbitrarily low "energy levels" for the perturbations. Thus the system is unstable and the 
instability is of a catastrophic nature since the increment \ fl\ has no upper bound. 

B. The case h > 0, Ui > 0; r— >Oatii— >oo,a scalar- type central singularity where 
the "scalar field" fj, is infinite. At rt — > 00 

^^0, V{x)^-^{l + o{l)). (65) 

The boundary condition at x — > is 

|y|/(-\/^lnx) < 00 (66) 

while the solution has the following asymptotic form at a; — > 0: 

y^V^iC^ + Celnx). (67) 

Again both the potential and the boundary conditions are d-independent. Considerations 
similar to those in item A lead to the same conclusion, namely, that the system is catas- 
trophically unstable. 

C. For the {d + 3)-dimensional section of the D-dimensional metric there is a range of 
integration constant values when the geometry has the form of a "space pocket" (see item 
C in Section 3). On the contrary, in our (0?+ 3) -dimensional metric (42), (43), containing an 
additional conformal factor, such a possibility is not realized. Indeed, using the expression 
(43) for (3 — \nr at the asymptotic u — > 00 and the relation (26) among the constants, one 
can verify that r can tend either to zero, or to a nonzero finite value corresponding to a 
black hole. 

D. Let us now consider perturbations of multidimensional black holes described by (20)- 
(24) under the conditions (31). The effective potential (59) may be written in an explicit 



15 



form using the transformation (28) in terms of the R coordinate: 



V{x) = Vi{x) + V2{x), 




d+1 R-2k 



-{{R + pN_)[R{2k + p - pN_) + 2kpN_] 



+{R - 2k) [Rp{l -N_) + ^{R + pN_f]}, 



2pN_{2k + pN_)R{R-2k) 



(68) 



where A^. = (A^ - d)/[N{d + 1)]. 

As follows from (68), Vi > and V2 > when R > 2k. The boundary conditions (54), 
(55) correspond to those conventional in quantum mechanics and thus the positiveness of 
V{x) means that solutions to (57) with < are absent. Consequently, multidimensional 
BHs are stable under monopole perturbations. Other types of multidimensional spherically 
symmetric solutions are strongly unstable. 

This generalizes the conclusions of Refs. [6,18] for the case of arbitrary d. 



We have seen that multidimensional black holes are stable under monopole perturbations, 
while other types of (electro)vacuum spherically symmetric solutions are strongly unsta- 
ble. This property distinguishes the black hole (BH) solutions from all the other possible 
solutions to the multidimensional field equations and supports the view that even if space, 
time and gravity are described by some multidimensional model, realistic collapse of iso- 
lated bodies should lead to black hole formation, as it is conventionally asserted in general 
relativity. In particular, this favours the models with a great number of primordial BHs 
present in the Early Universe, at epochs when there was no crucial distinction between the 
physical and internal dimensions. 
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